The stochastic entry of enveloped viruses: Fusion vs. endocytosis 
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Viral infection requires the binding of receptors on the target cell membrane to glycoproteins, or 
"spikes," on the viral membrane. The initial entry is usually classified as fusogenic or endocytotic. 
However, binding of viral spikes to cell surface receptors not only initiates the viral adhesion and the 
wrapping process necessary for internalization, but can simultaneously initiate direct fusion with the 
cell membrane. Both fusion and internalization have been observed to be viable pathways for many 
viruses. We develop a stochastic model for viral entry that incorporates a competition between 
receptor mediated fusion and endocytosis. The relative probabilities of fusion and endocytosis of a 
virus particle initially nonspecifically adsorbed on the host cell membrane are computed as functions 
of receptor concentration, binding strength, and number of spikes. We find different parameter 
regimes where the entry pathway probabilities can be analytically expressed. Experimental tests of 
our mechanistic hypotheses are proposed and discussed. 



I. INTRODUCTION 

Viral entry mechanisms are typically classified as either 
endocytotic, or as fusogenic In the latter, the virus 
membrane, after association with the surface of the host 
cell, fuses and becomes contiguous with the cell membrane. 
This process is mediated by the binding of cell surface recep- 
tors to glycoprotein spikes on the viral membrane surface, 
forming fusion competent complexes spanning the viral and 
cell membranes. In endocytosis, the host cell first inter- 
nalizes the virus particle, wrapping it in a vesicle before 
acidification-induced fusion with the endosomal membrane 
can occur. Wrapping can occur only after cell surface recep- 
tors, that also act as attachment factors, bind to the viral 
spikes. Experimentally, both fusion with the cell membrane 
and internalization can be observed and distinguished using 
microscopy 0, Q • 

Many viruses, such as influenza and hepatitis B, use endo- 
cytosis as their primary mode of entry [4j, |5J . The surface of 
an influenza virus is coated with ~ 400 hemagglutinin (HA) 
protein spikes @. The HA adheres to sialic acid-containing 
glycoproteins and lipids on the cell surface leading to wrap- 
ping of the virus particle. Particle wrapping may also be me- 
diated by the recruitment of pit-forming clathrin/caveolin 
compounds 0]. Cell membrane pinch-off, leading to inter- 
nalization of the virus, usually requires additional enzymes 
such as dynamin and endophilin Q. Endosomal acidifica- 
tion oligomerizes the HA, priming them to fuse with the 
endosomal membrane. Direct fusion of the influenza virus 
with the host cell membrane is precluded since HA is acti- 
vated only in the acidic endosomal environment. However, 
low pH conditions have also been shown to induce the direct 
fusion of influenza virus with certain cells [Iol |. 

Recent experiments on the avian leukosis retrovirus 
(ALV) have provided evidence both for a pH-dependent di- 
rect fusion mechanism [ill . [l2| , and an endocytotic pathway 
[l3l |. Moreover, the entry pathway of some viruses such as 
Scmliki Forest Virus (SFV) can be shifted from endocyto- 
sis to fusion by acid treatment, but only in certain host cell 
types [3]. In this case, low pH triggering of receptor-primed 
envelope glycoproteins can initiate fusion before the virus 



can be wrapped and endocytosed. Vaccinia and HIV (typi- 
cally infecting cells via direct fusion) have also been shown 
to exploit both entry mechanisms @, HH, Oil Gj} ■ For exam- 
ple, fusion-independent mechanisms of HIV-1 capture and 
internalization in mature dendritic cells, mediated by DC- 
SIGN [lH, can be a significant mode of HIV transmission 
through dendritic cells and lymphatic tissue (l9l |. Capture 
by DC-SIGN and CLEC-2 adhesion molecules also internal- 
izes HIV in platelets [2(|. Thus, depending upon physical 
conditions and cell type, both entry pathways arc potentially 
accessible to certain viruses. The choice seems to depend on 
the type of receptors the viruses engages, whether they are 
receptors/coreceptors that induce fusion (perhaps triggered 
by low pH), or simply attachment factors such as sialic acid- 
rich glycoproteins that do not induce fusion. In this latter 
case, complete wrapping before an irreversible fusion event 
is more likely to occur, and internalization is favored. 



It is not surprising that subtle changes in the interactions 
between viral membrane proteins and cell receptors dramat- 
ically affect the infectivity of a virus, as recently demon- 
strated for the 1918 influenza virus [21(. In this paper, we 
model virus-receptor kinetics and propose a mechanism con- 
sistent with the above experimental observations, and that 
describes viral entry by incorporating both fusion and en- 
docytosis entry pathways in a probabilistic manner. In the 
next section, we develop a stochastic one-species receptor 
model for the binding of receptors necessary to start the 
virus wrapping process. These receptors, upon binding, can 
induce membrane fusion at each receptor-spike complex. In 
the Results, wc find parameter regimes in which each of 
the entry mechanisms dominate. Explicit expressions for 
the entry pathway probabilities, as functions of the relevant 
kinetic rates, are given in the Appendix. In the Discus- 
sion and Summary, we explore the connection between our 
parameters and experimentally controllable physical condi- 
tions such as rcccptor/coreceptor density, spike density, and 
cell membrane rigidity. Experimental tests are proposed and 
extensions of our analysis to more realistically incorporate 
biological features are discussed. 
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FIG. 1: A schematic of viral entry pathways, (a) A nonspecifically adsorbed virus particle can desorb with rate kd, or it can (b) 
recruit and specifically bind a receptor. The receptor can immediately initiate membrane fusion with rate kf as shown in (c), or, it 
can recruit additional receptor molecules, inducing wrapping of the virus particle. From partially wrapped states (d), the virus can 
at any stage undergo membrane fusion (e), or, it can completely wrap and internalize the virus particle ((f) and (g)). 



II. SINGLE RECEPTOR KINETIC MODEL 



The basic features of our proposed mechanism are shown 
in Fig. [TJ A virus particle initially nonspecifically adsorbed 
on the cell membrane (without any bound receptors) can 
spontaneously dissociate with rate kd- Alternatively, mo- 
bile receptors in the cell membrane can bind specifically to 
the glycoprotein spikes (assumed here to be uniformly dis- 
tributed) on the viral surface. Successive addition of re- 
ceptors to the viral ligands, when n are already attached, 
occurs with rate p n . Thus, the binding of the first receptor 
occurs with rate pq. Similarly, desorption of the n th recep- 
tor occurs at rate q n . We consider the adsorption of a single 
effective receptor or attachment factor to a spike, lumping 
together the effects of multiple reccptor/coreceptor types. 
This approximation is valid when, for example, coreceptor 
binding is highly cooperative such as suggested in the HIV 
infection process where CD4 binding to the gpl 20 protein 
spike induces rapid CCR5 coreceptor binding [22 |. 

Receptors not only adhere the cell membrane to the vi- 
ral membrane, but can also initiate local membrane fusion 
at each receptor-spike complex with rate kf. Fusion can 
occur at any time during the receptor recruitment process 
and is more likely to occur per unit time with more bound 
receptors. Rcceptor-spikc complexes that are unable to ini- 
tiate fusion arc described by a vanishing fusion rate kf — > 0. 
However, if receptors have high fusogenicity kf, the virus 
might fuse only after a single receptor has attached. Only 
if the system reaches a fully wrapped state with N bound 



receptors (Fig. [Tf), before any fusion event occurs, can 
pinch-off and endocytosis occur with rate k e . The num- 
ber of spikes A^ is ty pica lly large, varying among viruses 
such as HIV (N ~ 15) (H, SIV (N ~ 70)01, and Influenza 
(N ~ 400)0]. The path to endocytosis is thus a race be- 
tween fusion and complete wrapping of a relatively large 
number A" of viral spikes. 
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FIG. 2: The stochastic process representing the competition be- 
tween membrane fusion and endocytosis. The states n corre- 
spond to the number of receptor-spike complexes formed, while 
N is the total number of spikes on the virus membrane. Each 
receptor-spike complex can initiate membrane fusion with rate 
kf. As more receptors are bound, the total rate of fusion in- 
creases linearly. The irreversible pinch-off and endocytosis rate 
is denoted k e . 

States in the model are labeled by the index n (Fig. 
representing the number of formed receptor-spike com- 
plexes. Starting from a nonspecifically adsorbed virus parti- 
cle denoted by state n = 0, the system progresses along the 
chain with the appropriate transition rates corresponding 
to attachment and detachment of cell membrane receptors. 
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The probability P n (t) of having n bound receptors (or at- 
tachment factors) at time t obeys the master equation 

Pn(t) = ~{nkf +p n + q n )P„ +p n -lPn-l 

+q n+1 P n+ i, 1 < n < N - 1, 

(1) 

Po(t) = -{kd+Po)Po + qiPi, 

P N (t) = -(Nk f + k e + q N )P N + p N -iP N --i_. 

Although we describe the general viral entry process in 
terms of recruitment and binding of a single type of re- 
ceptor, our model encompasses processes involving clathrin 
or calveolin aggregation and pit formation, typically lead- 
ing to endocytosis. All of these biologically distinct, but 
physically similar mechanisms can be analyzed by appropri- 
ately interpreting the rates. For example, the nucleation of 
a clathrin-coated pit can be modeled by effective binding 
and unbinding rates p n and q n that describe the rates of 
clathrin addition and removal from the pit. Since individual 
clathrin molecules are not known to induce membrane fu- 
sion, the fusion rate kf = 0, and only endocytosis (or virus 
dissociation from the cell surface) would occur. The fusion 
rate is also negligible if the receptor is a simple attachment 
factor that adheres the membranes, but does not facilitate 
fusion. 

The receptor binding and unbinding rates, p n and g n , are 
related to the cell surface receptor density and the receptor- 
spike binding strength, respectively. The Markov process 
shown in Fig. [2] implicitly assumes that the receptor re- 
cruitment is not diffusion-limited - the rate of addition of 
successive receptors is independent of the history of previous 
receptor bindings. 

Approximate forms for p n , q n can be physically motivated 
by considering the number of ways additional receptors can 
bind or unbind, given that there already exist n receptor- 
spike complexes. As the membrane progressively wraps 
around the virus particle, the rate of addition of the next re- 
ceptor is proportional to the number n p of unattached spikes 
bordering the virus-cell membrane contact line L, as shown 
in Fig. O 

For a large number TV of uniformly distributed spikes 
on a virus particle of radius R, the contact area A c « 
^R 2 n/N w 2ttR 2 {\ - cos6»„) where 0„ = cos" 1 [1 - 2n/N] 
is the angle subtended by the contact line when n receptors 
are attached. Since the area per spike is a s ~ AttR 2 /N , the 
number of spikes near the contact perimeter L = 2ttR sin n 
is found from n p f=a L/ ^/aj. Upon using the explicit form for 
9 n , we find the number n p of periphery spikes as a function of 

n > 1 receptor-bound spikes, n p ~ y/nN \J 1 — (1 — 2n/N) 2 . 
At a stage where n > 1 receptors have bound, approximately 
n p spikes are accessible for additional binding of receptors. 
Similarly, there are approximately n p receptor-spike com- 
plexes available for dissociation. Combinatorically, the re- 
ceptor binding and unbinding rates take the form [24| 
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FIG. 3: Schematic of a partially wrapped virus particle. The 
unbound spikes above the contact region are represented by open 
circles, while the receptor-bound spikes in the contact region are 
represented by the red-filled circles. The number n p of spikes 
or spike-receptor complexes near the contact perimeter used to 
compute p n or q n via Eqs. [5] are shown as black dots. 



Pn~Po—i= q n ^ qi —= (2) 

where 1 < n < N — 1, po is the intrinsic rate of binding the 
first receptor, and q\ = q^ is the dissociation rate of an in- 
dividual receptor-spike complex. Since q\ is a spontaneous 
receptor-spike dissociation rate, it is independent of receptor 
density and spike number. If the receptor-spike binding en- 
ergy is at least a few fcgT, we also expect q n to be relatively 
insensitive to the cell membrane bending rigidity. 

A number of physical attributes and biological interme- 
diates can be incorporated into the rate parameters to ad- 
dress more complicated microscopic processes. For example, 
if thermal fluctuations are rate-limiting, there would be an 
additional factor in p n reflecting the probability per unit 
time that a patch of membrane fluctuates to within a dis- 
tance of the virus surface spike that allows receptor-spike 
binding. The dynamics of this process depends on the cell 
membrane tension and rigidity, the typical spike spacing, 
and potentially the viscosity of the extracellular environ- 
ment. For stiff membranes under tension, the wrapping of 
spherical par ticles encounters an energy barrier near half- 
wrapping [25| , which can be incorporated into the dynamics 
by assuming an additional factor in p n that has a minimum 
near n ~ N/2. Thus, the energy barrier associated with 
membrane bending will tend to flatten the n-dcpcndcnce of 

Pn- 

Particular aspects of the entry process can also influence 
estimates of the other rates. If the viral spikes are mobile 
and can aggregate to the initial focal point of adhesion, the 
target cell membrane is not able to fully wrap and endocy- 
tosis is prevented. The overall kinetic scheme remains un- 
changed except with k e « 0. Finally, the recruitment of sec- 
ondary coreceptors that occurs in, e.g. HIV fusion, can also 
be developed within our current framework and is discussed 
in the Discussion and Summary section. Although some of 
the physical details described above influence the specific 
values of p n , we will show that for large iV, the qualitative 
behavior of our model can be summarized by distinct pa- 
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FIG. 4: The currents through each pathway for TV = 20 spikes. 
The dashed black curve is the current for desorption, while the 
red and green curves are the currents for fusion and endocytosis, 
respectively. Time is measured in units of k^ 1 , and all rates are 
normalized with respect to kd- The parameters used in all plots 
are po = 1, qi = 0.3, and k e = 0.3. (a) The currents for constant 
Pn = 1, q n = 0.3 and fusion rate kf = 0.001. (b) The same 
parameters except that Eqs. [2] are used for the rates p„,q n . (c) 
The currents with p n , q n as in (b), except that the fusion rate of 
each spike-receptor complex is increased to kf = 0.01. 



rameter regimes, somewhat insensitive to the precise values 
of p„,q n - 
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FIG. 5: Numerical solutions of entry probabilities Qi (all rates 
normalized by kd)- (a) The entry probabilities as a function 
of po- Endocytosis arises only for larger po > qi, after the 
fusion probability becomes significant. Parameters used are 
TV = 20, kf = 0.003, fee = 0.3 (all rates are normalized by kd). 
(b) The probabilities of dissociation (dashed), fusion (red), and 
endocytosis (green) as functions of the individual receptor-spike 
fusion rate kf. Here, qi = 0.3. The thin dashed red curve cor- 
responds to a faster receptor detachment rate (qi = 2) which 
prevents endocytosis. 



III. RESULTS 

We solved Eq. [T] for the probabilities P„ (t) given an ini- 
tial condition P n (t = 0) = S n ,o- Using these probabili- 
ties, we find the probability currents through the dissoci- 
ation, fusion, and endocytosis pathways, Jd{t) = kdPo(t), 

Jf(t) = ^/Eri^W. and J e(t) = k e P N (t), respectively. 
Figure 3] shows the current for each pathway as a function 
of time (in units of k^ 1 ). Note that detachment, fusion, 
and endocytosis arise sequentially in time. Upon comparing 
FigHJa) and (b), we see that changing p n , q n from constant 
values to the forms in Eqs. [2] shifts the currents through 
the fusion and endocytotic pathways to earlier times. This 



speed-up is simply a consequence of the larger hopping rates 
Pn,qn, especially for n ss N/2. Nonetheless, the specific 
form of p n ,q n , provided they are slowly varying in n, only 
quantitatively affect the timing of the onset of the currents. 
The dramatic variations in the infection pathway taken come 
with changes in kf. When the fusion rate kf is increased, 
endocytosis is suppressed in favor of fusion as shown by Figs. 
Hb) and He). 

Upon time-integrating the currents, we find the total 
probabilities Qi = J Q Ji{t)At, (i = d, e, /) for each pathway. 
Note that from probability conservation, Qd + Q e + Q f = 1. 
In Figure [51 we use the binding and unbinding rates given 
by Eqs. [2] to numerically compute the entry pathway prob- 
abilities Qi. Fig. [5ja) shows the pathway probabilities as a 
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function of the intrinsic receptor binding rate po/qi- This 
ratio is a measure of the density-dependent free energy AG 
of the spike-receptor binding: po/qi ~ e AG / kT [24)]. In the 
simplest limit of extremely low receptor density, (po — + 0), 
Qd ~ 1 - O(p ),Qf ~ O(kfpo), and Q e ~ 0, and only 
dissociation can occur. As po is increased, the probability 
of fusion increases at the expense of desorption. Endocy- 
tosis remains negligible as long the states that occur with 
any appreciable probability are those with small n. Only as 
Po ^> qi does the probability of endocytosis become appre- 
ciable and approach the asymptotic expression given by Eq. 
[3] in the Appendix. 

Figure E^b) shows the total probabilities Qd of dissoci- 
ation, Qf of fusion, and Q e of endocytosis as functions 
of the fusion rate kf. For large detachment rates (e.g. 
qi = 2 > po = 1), Q e ~ 0, and fusion can occur only at large 
kf, as shown by the thin dotted curve. For the parameters 
used, TV = 20, po = 1, <Zi = k e = 0.3 (normalized by kd), the 
transition from a predominantly endocytotic pathway to a 
predominantly fusion pathway occurs for 10 -3 < kf < 10 -2 . 
When kf 3> 10~ 2 , the sum of fusion probabilities over all 
intermediate states is appreciable, preventing endocytosis. 
Therefore, only for a particularly small fusion rate kf, and 
nonnegligible endocytosis rate k e is internalization possible. 
We show in the Appendix that generally, if p n > q n , the 
effective drift of the stochastic system towards the wrapped 
state renders the partitioning between fusion and endocy- 
tosis controlled by the fusion rate kf. If this is the case, 
the transition from endocytosis to fusion occurs at approxi- 
mately kf ~ (p n + q n )/N 2 , as is confirmed by the numerical 
solution for TV = 20 shown in Fig. EJb). 

Although we have chosen TV = 20 as a repreentative 
spike number, different viruses and their variants can have 
a widely varying number of active spikes. In Fig. [5] we show 
how the entry pathway probabilities depend on the number 
TV of active viral spikes. As TV is increased, the probabil- 
ity for fusion increases at the expense of endocytosis. For 
TV — > 00, and a nonzero kf, Q e — * since fusion will likely 
occur during the infinitely long wrapping processes. For 
small TV, the plotted probabilities must be interpreted with 
an TV-dependent k e . Suppose TV < 10. Even if all spikes are 
receptor-bound, the membrane has an appreciable distance 
to bend and before full wrapping and endocytosis can oc- 
cur. Effectively, the fusion rate k e starts to decrease if TV 
gets small such that the rate of membrane fluctuations over 
a typical interspikc distance decreases. 

Provided p n — q n 1/TV, asymptotic analysis of the 
solutions reveal qualitatively different behaviors depending 
upon how the fusion rate kf compares with 1/TV. If the 
receptor-spike complex is highly fusion competent such that 
kf/(pn + qn) 3> 1/TV, the probability of reaching a com- 
pletely wrapped (n f=a TV) state is exponentially small and 
endocytosis cannot occur. Here, the virus pathway is nearly 
entirely partitioned between dissociation and fusion, as in- 
dicated by the asymptotic expressions for Qi given by Eqs. 
2] in the Appendix. 

If 1/TV 2 <C kf/(p n + q n ) <C 1/TV , the receptor-spike com- 
plex has intermediate fusogenicity. In this case, the time- 
integrated probability J Q P n (t)dt used to construct Qi re- 
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FIG. 6: The pathway dependence on receptor associa- 
tion/dissociation rates and the number TV of virus spikes. The 
number of spikes controls which regime Eq. [4] or Eq. [5] is valid. 
Large TV enhances fusion almost entirely at the expense of endo- 
cytosis. 

mains small for n ~ TV and the probability of endocytosis 
is still exponentially small, despite the smaller kf. In this 
regime, we find in the Appendix (Eq. [5]), an expression for 
Qd- Only when kf/(p n + q n ) <C 1/TV 2 is virus internaliza- 
tion appreciable. Expressions for Qf and Q e in this very 
weak fusogenicity limit are also displayed in Eq. [5] of the 
Appendix. The expression we find agree with the exact nu- 
merical evaluation of Qd.ej from solving Eq. [T] 

IV. DISCUSSION AND SUMMARY 

Fusion can be directly distinguished from endocytosis by 
imaging fluorescent markers loaded into viruses or model 
vesicles. Upon fusion, one expects to see an immediate re- 
lease of marker into the periphery of the target cell. Simi- 
larly, singlc-liposomc fluorescence imaging experiments that 
label and detect vesicle lipids and as they mix with a sup- 
ported bilayer upon fusion (26j , can be used as in vitro model 
systems for virus fusion. In such experiments, the model 
parameters p n ,q n ,kf,k e ,kd can be tuned by controlling cer- 
tain physical chemical properties, enabling one to dissect the 
mechanism of viral entry. The entry pathways delineated 
by the different parameter regimes described in the previ- 
ous section, and by the asymptotic formulae given in the 
Appendix, provide a framework for analyzing and designing 
viral entry experiments. 

The receptor density plays the first critical role via the 
binding rate parameter p n . For low receptor densities, and 
proportionately lower p n (but fixed spontaneous detachment 
rate q n ), the virus particle can only dissociate or fuse. Al- 
though lowering receptor concentration decreases the overall 
entry probability, it can increase Qf/Q e , the fusion probabil- 
ity relative to endocytosis probability (see Fig. 02a)). Only 
for p n — q n > 0(1/ N) can the receptor spike complex fusion 
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rate k f become important in determining whether fusion or 
cndocytosis occurs. In order for endocytosis to occur, the 
fusion rates must be small such that kf/(p n + q n ) <C 1/N 2 . 
The rate p n can also be substantially decreased by increas- 
ing the target membrane surface tension, thereby suppress- 
ing the thermal fluctuations of the membrane required to 
bring cell receptors and viral spikes into proximity. 



The rapid drop-off in endocytosis predicted as the fusion 
rate is increased from kf/(p n + q n ) -C 1/TV 2 to kf/(p n + 
?n) ^ 1/^ 2 j especially for large TV, shows that tuning 
physical conditions (such as pH or temperature) that af- 
fect the fusogenicity of receptor-spike complexes, kf, can 
have a large effect on the viral entry pathway. Recent ex- 
periments by Melikyan et al. [27[, Henderson and Hope 
[281 ] . and others have shown a rate-limiting intermediate in 
the HIV fusion process that can be arrested by lowering 
temperature. Since CD4 binding was not the rate limiting 
step, lowering the temperature decreases kf to a degree pre- 
sumably much less that (p n + q n )/N 2 , preventing fusion. 
If these systems have the necessary endocytotic machinery 
and support pinch-off, lowering temperature and arresting 
the receptor-spike fusion complex while retaining the adhe- 
sive wrapping of receptor-spike binding would enhance the 
endocytotic pathway. The effective rate kf can also be low- 
ered by cross-linking (with e.g., defensins) membrane glyco- 
proteins, rendering their complexes with viral spikes fusion 
incompetent [291 ] - However, if the cross-linked glycoproteins 
retain their attraction for the viral surface, the probability 
of wrapping and internalization would increase. If an inde- 
pendent measurement or estimate of p n +q n is available, the 
dependence of kf on temperature, pH, and chemical modi- 
fication can be probed. 



Finally, pathways to fusion and endocytosis can diverge 
for systems that require both primary receptors and sec- 
ondary receptors (coreceptors). If two species are required, 
one for adhering cell and viral membranes, and another 
to induce fusion, endocytosis will be favored, all else be- 
ing equal. In this case, the initial receptor binding only 
causes the cell membrane to wrap around the virus parti- 
cle. An additional corcccptor must diffuse and bind to the 
spike-receptor adhesion complex to induce fusion. A highly 
cooperative receptor-coreceptor interaction, such as in HIV 
fusion involving CD4 and CCR5/CXCR4, is still modeled by 
Eq. [TJ but with the binding rates p n interpreted as an effec- 
tive binding rate for both classes of receptor, proportional 
to the product of their surface concentrations. However, 
if the coreceptor density and/or mobility is limiting [l7| . 
the binding of receptors occur first, with coreceptor priming 
and formation of a fusion competent receptor-spike complex 
occurring slowly. This allows time for receptor (adhesion 
molecule) mediated membrane wrapping of the entire virus, 
enhancing the likelihood of endocytosis. Thus, by main- 
taining a high adhesion receptor density, and lowering the 
fusion-enabling coreceptor density, one enhances the endo- 
cytotic pathway. 



V. APPENDIX 

Consider the Master Equation [1] in the form P(t) = 
— [Mo + M/] P, where Mo is the conserved random walk 
transition matrix involving only kd,p n ,qni an d k e , and 
Mj = fc/diag(n) is decay term arising from fusion. Large 
TV expressions for the entry pathway probabilities Qi can be 
obtained in different limits. 

In the limit where the receptor binding is irreversible 
{Pn/q-n — * oo), Eq. Q]can be solved exactly: 



Qd = 
Q f = 

Qe = 



Pa + k d 
1 

Po + k d 



kePN-l TJ 



Pm-1 



fee + Nkf £± p m + TJlkf 



(3) 



PN-1 TT Pm-l 



n 



fe e + Nkf po + kd p m + mkf 



If q„ =0, the probability of dissociation is fixed by k d and 
Po, and the remaining current is partitioned between fusion 
and endocytosis. For a given po,kf,k e , this limit gives the 
highest probability of the maximally-wrapped state and the 
highest endocytosis probability. 

When p n — q n ^> 1/TV, there are two possible limits corre- 
sponding to high receptor-spike fusion rates, kf / (p n +q n ) ^> 
1/TV, and intermediate fusion rates 1/TV 3> kf/(p n + q n ) ^> 
1/N 2 . For k f /(p n + q n ) > 1/TV, f °° P„(t)dt is nonnegligible 
only for n ~ and J Q °° Pjy(t)dt is exponentially small. A 
small n local analysis of Eq. Q] yields 



Qd 
Q f 



kd + po- poqi/ (pi + qi) 
Po 

kd +Po + k d q\/p\ 



and 



(4) 



0. 



for J- — > 1/TV, 

p n + q n 



explicitly showing the absence of endocytosis. The values 
for Qd and Q / corresponding to the parameter values used 
in Fig. [6] are Q d ~ 0.56 and Qf as 0.435, which agree well 
with the large TV limit of the numerical soution. 

In an intermediate fusogenicity regime, 1/TV 3> kf/(p n + 
q n ) ^ 1/N 2 , J °°P„(i)d< is nearly constant for n <C TV. 
However, J °°P Ar (i)dt remains small and endocytosis, al- 
though still unlikely, occurs with slightly higher probabil- 
ity than in the high fusogenicity limit. In this intermediate 
limit, one might first attempt perturbation theory about 
kf = (since its largest element of M/, Nkf -^i p n + q n , 
is much small than the typical elements of Mo). However, 
as 1/TV decreases, so do the spacings between eigenvalues of 
the zero- fusion transition matrix Mo. Heuristically, for per- 
turbation theory to be accurate for slowly varying pn^n, 
the largest of the diagonal correction terms, Nkf, must be 
smaller than (p n + q n )/N, the typical spacing between eigen- 
values. First order perturbation for all P n (t) is accurate 
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only if kf/(p n + q n ) *C 1/N 2 , as is explicitly shown by 
Eqs. [3] Upon expanding Q/ and Q e (from Eqs. [3]) in 
powers of kf, one finds 0(N 2 ) corrections terms. Thus, the 
expansion is only accurate if kf/(p n + q n ) <C 1/N 2 . En- 
docytosis is preempted by fusion or dissociation only when 
perturbation theory about a nonnegligible Q e fails (when 
kf/{Pn + In) it 1/N 2 ). Perturbation results for the few- 



receptor states n w important for the dissociation proba- 
bility Qd, are still valid as long as kf /(p n +q n ) -C 1/N. How- 
ever, the perturbation results that include the flux through 
larger (n w AT) states are accurate only if the receptor-spike 
complexes are extremely inefficient at initiating membrane 
fusion, and kf/(p n + q n ) <C 1/N 2 . Thus, for small enough 
kf/{p„ +q n ), we find 



Qd 
Q f 

Qe : 



k d + k d Y,n=i U^=i(n/Pi) + (k d q N /k e ) njl^fa/Pi) 

kd + PO + k d J2n=l lir=l fe/P») + (kdqw/ke) YliLl'ili/Pi) 



kf/(p n + q n )<N- 1 



k e+P ^k d k e [i + Ef^ 1 m = i (*)] + ^nS 1 



fe //(Pn + 9n) < AT 2 where 
, k f /(p n + q n )<N- 2 , 



(5) 



independent of kf. Equations[3][S]givc estimates for the entry 
probabilities Qd.ej in the different parameter regimes. For 
small N such that fe//(p n + q n ) < A^ -2 , Q e w 0.411, which 
agrees well with the limit shown in Fig. [6] 
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